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Abstract. We report the first result of a supernova search program designed to measure the evolution of the 
supernova rate with redshift. To make the comparison with local rates more significant we copied, as much as 
possible, the same computation recepies as for the measurements of local rates. Moreover, we exploited the mul- 
ticolor images and the photometric redshift technique to characterize the galaxy sample and accurately estimate 
the detection efficiency. 

Combining our data with the recently published meaurements of the SN la rate at different redshifts, we derived 
the first, direct measurement of the core collapse supernova rate at z = 0.26 as r cc = 1.45^0 45 h 2 SNu [/i=Ho/75]. 
This is a factor three (± 50%) larger than the local estimate. The increase for a look back time of "only" 2.8 
Gyr is more rapid than predicted by most of the published models of the SN rate evolution. Core-collapse SN 
rates measure the death rate of massive star and, because of the short time scale of evolution, can be translated 
in a measurement of the ongoing SFR. Assuming a Salpeter IMF and the standard scenario for core-collapse 
progenitors we derived to an estimate of the star formation rate at redshift 3.1_ J' x 10 -2 /i 3 M0yr _1 Mpc -3 
which compare very well with a recent estimate based on the measurement of the Ha luminosity density at the 
same redshift. 

Key words. supernovae:general - star:formation - galaxy:evolution - galaxy:stellar content 



1. Introduction 

Supernova rates represent a link between the evolution of 
individual stars and that of stellar systems. In particular, 
the rate of type Il + Ib/c SNe measures the death rate 
for core-collapse (CC) of young, massive stars and, be- 
cause of the short time-scale of the progenitor evolution 
jHeger et al. 1200*31) , directly reflects the on-going star for- 

Send offprint requests to: E. Cappellaro, e-mail: 
cappellaro@na . astro . it 

* Based on observations collected at the European Southern 
Observatory, Chile (ESO Programmes 62.H-0833, 63.H-0322, 
64.H-0390, 67.D-0422, 68.D-0273, 69.D-0453). 



mation rate (SFR) in a given environment. On the con- 
trary, the rate of type la S Ne, which result from long-lived, 
low mass binary systems llBranch et al. lfl995|) reflects the 
long-term star formation history. Owing to the very high 
intrinsic luminosity, SNe can be seen at very large dis- 
tances and hence measurements of SN rate evolution with 
redshift can be used to trace the history of the SFR with 
cosmic age. 

Studies of the evolution of SN rates with resdhift are 
now strongly unbalanced towards theory. In recent years, 
many authors have published predictions of the SN rate as 
a function of redshift based on the SN progenitor scenar- 
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ios a nd the modeling of the cosmic star formation his- 
tory llMadau. Delia Valle fe Panaeial Il998t ISadat et al. I 
1998tlDahlen fc Fransson ll999HYungelson fc Liviol2000l: 
Ko bayashi. Tsuiimoto. fe Nomotd l2000t ISullivan et al. I 
2000t ICalura fc Matteucci Il2003j) . Despite the strong in- 
terest, observational estimates of the SN rate at high red- 
shift are still very scanty. The few published measurements 
are based on SN sea rches aimed at using type la SNe as 
cosmological probes l)Pain et al. Ill99fil l2002t iTonrv et al. I 
120031) and, as a consequence, they are strongly biased to- 
wards type la. To date, there is no direct measurement of 
the evolution of core-collapse SN rates. 

With the goal to fill this gap, we initiated a long term 
project aimed to measure the rate evolution with redshift 
for all SN types. In this paper we report the first results 
of this effort, namely an estimate of the core-collapse SN 
rate at redshift z ~ 0.3. 

To reduce the systematics in the comparison with the 
local rate, our strategy we followed as c losely as possible 
the s a me approach used in that context JCapoellaro et al.1 
Il997t ICappellaro. Evans, fc Turattol Il9991. The process 
consists of four steps: i) SN candidate detection and, when 
possible, spectroscopic classification (Sec. 121; it) charac- 
terization of the galaxy sample through galaxy photomet- 
ric rcdshifts; estimate of distances and calibration of the 
absolute luminosities for each galaxy (Sec.EJ ; in) evalua- 
tion of the SN detection efficiency and thus of the effective 
surveillance time for each galaxy of the sample (Sec. HJ, 
iv) estimate of the SN rate per unit luminosity. To round 
up the paper, we discuss the main sources of uncertainty 
(Sec. 6) and discuss the implication of our results (Sec. 7). 

2. The supernova search 

The building block of a SN search is the detection of vari- 
able sources through the comparison of images of selected 
sky fields obtained at different epochs. In general, the tem- 
poral sampling of the observations is tuned to the specific 
goal one wants to achieve. For the use of type la as cos- 
mological distance indicators it is crucial to catch SNe 
as early as possible and hence, accounting for the typical 
rise time, the observations have to be spaced by 2-3 weeks. 
Instead, to maximize the event statistics, the time elapsed 
between exposures should be equal (or longer) than the 
time a typical event remains brighter than the search de- 
tection limit. The latter, of course, depends on the target 
distance (or redshift) and the SN type. Independently 
on the temporal sampling, to ensure that all SNe 
are detected, the time elapsed between the first 
and the last observation of a given field have to be 
longer that time for a significant luminosity evo- 
lution for all SN type. The latter can be as long 
as 3-4 months for SN IIP and even longer for type 
Iln. 

For our search we selected 21 fields, evenly distributed 
in right ascension, which have been monitored for about 2 
years with an average sampling of one observation every 
three months. However, the results reported in this paper 



are based on observations of 5 fields only, i.e. those with 
the best temporal and filter coverage. The full sample will 
be presented in a future paper (Riello et al. 2004, in prepa- 
ration). The observing log is shown in Tab. ^ For 
each field we give the center coordinates and, for 
the epochs when observations are available, we list 
the seeing. 

A typical observing run was splitted into two parts: 
the search and the follow-up observation of candidates. 

For the search, two consecutive nights were devoted at 
the ESO/MPI 2.2m telescope at ESO, La Silla (Chile). 
The telescope was equipped with the Wide Field Imager 
(WFI) and a mosaic of 2 x 4 CCD detectors of 2048 x 4096 
pixels which image a sky area of ~ 0.25 deg 2 with an 
excellent spatial resolution of 0.238 arcsec/pix. 

When possible, the first observing night was dedi- 
cated to obtain deep V band exposures for candidate de- 
tection while in the second night the same fields were ob- 
served through a different filter, B or R, with the purpose 
of collecting color information both for the candidates and 
the galaxies. Unfortunately, due to a number of tech- 
nical, meteorological and scheduling constraints in 
many cases we could not maintain this observing 
strategy as can be seen from Tab. ^ This implies 
that only in a few cases we could derive the candi- 
date color. For homogeneity, in the following sta- 
tistical computation we considered only the candi- 
dates detected in the V band exposures. 

In order to get rid of detector cosmetic defects, cosmic 
rays, satellite tracks and fast moving objects, for each field 
we obtained three 900 s exposures dithered by 5-10 arcsec. 

Follow-up observations were scheduled about one week 
after the search at the VLT+FORS1/2 at ESO Paranal for 
the spectroscopic classification of some of the candidates. 
The VLT was needed, as most of our SN candidates are 
in the magnitude range V ~ 22.5-23.5 mag. For a proper 
subtraction of the night sky emissions we selected grisms 
of moderate resolution, namely grism 300V and/or 3001 
(resolution ~ 10A FWHM), which allowed us to cover a 
quite wide wavelength rang e (4000 -r 11 000 A). Depending 
on the candidate magnitude, exposure times ranged from 
900s to 3 hours. Details on the reduction of the spectro- 
scopic observations and of the spectral analysis will be 
given elsewhere. Here we make use only of the spectral 
classification and redshift. 

The analysis of the search images began with the re- 
moval of the instrument signature and calibration for 
which we used IRAF 1 and MSCRED, a package specif - 
ically designed to handle mosaic images llValdeslll99l . 
Indeed, after bias subtraction and flat fielding, the indi- 
vidual dithered exposures were astrometrically calibrated 
and stacked in a single image. For photometric nights, 
observations of standard fields were used for absolute cal- 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 
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Table 1. Observing log of the SN search fields. For each field, when observations in a given band are available, we 
report the measure of the seeing in arcsec (FWHM on the stellar objects on the combined image) 
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ibration llLandolt 11 19921. Otherwise, the photometric zero 
point was established by comparison with a calibrated im- 
age of the same field. 

For each field we computed the difference between the 
image to be searched (target image) and a suitable archive 
frame (template image). Indeed, after accurate astromet- 
ric and photometric registrations, the most crucial step 
in this process is the matching of the point spread func- 
tion (PSF) of th e two image s. This was done using the 
ISIS2.1 package llAlard 11200(1 that, from the comparison 
of the same sources in the two images, computes a space- 
varying convolution kernel to degrade the image with the 
best seeing to match the other one. Taking into account 
that the best subtraction is obtained when the two images 
have similar PSFs, and that we want to preserve as much 
as possible the resolution of the target image, we had to 
maintain an archive with template images with different 
seeing. Populating the archive required a significant in- 
vestment of telescope time and this was the reason why 
the search became really efficient only after some time 
from the first observations. 

Variable sources leave residuals in the difference im- 
age which have been detected and l ogged into a cata- 
logue using the SExtractor program (Bert in & Arnouts 
1996) which has also the capability of separating stars 
from galaxies. Due to residuals of poorly removed bright 
stars or cosmic rays, the variable source catalogue contains 
many false detections, most of which are quickly elimi- 
nated by means of a custom made ranking program. This 
makes use of information from the difference image as well 
as from the target and template images and it has been 
tuned through extensive artificial stars experiments. The 
surviving candidates, typically a few tens per field, are all 



checked visually by a human expert. Among these, a few 
are still obviuos false detections which could not be prop- 
erly flagged by our software but are quickly eliminated 
by visual inspection. Among these are residuals of moving 
objects which are not completely masked by our dithering 
strategy. After that, we are left with true variable source 
among which we remove variable stars, objects with stel- 
lar profile present both in the target and template image, 
but with different magnitudes. At the end of this process, 
one is left with the fiduciary SN candidates, typically from 
none to a handful per field. 

Ideally, one would need spectroscopic confirmation for 
all the candidates. Unfortunately, because of the limited 
VLT time available, we could obtain spectroscopic obser- 
vations for ~ 20% of the detected SN candidates only. This 
is the main weakness of the work presented here. On the 
other hand, we could verify the reliability of our SN candi- 
date selection criteria: out of the 29 candidates for which 
we have obtained VLT spectra during the entire search 
project, 22 turned out to be SNe (45% type la and the 
other type II and Ib/c) and 7 variable AGNs. We stress 
that our approach to candidate selection was designed to 
avoid as much as possible any selection bias and in partic- 
ular, we do not exclude a priori nuclear candidates. Given 
that, the intrusion of AGNs is unavoidable. 

Even without spectroscopy, contamination by variable 
AGNs can be reduced by looking at the long term vari- 
ability history of the candidates. With this aim we kept a 
database, powered by MySQL 2 , which is used to search 

2 MySQL is an open source database without a license fee 
under the GNU General Public License (GPL). See the project 
homepage for further details, http://www.mysql.com/ 
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Table 2. SN candidates of the fields of Tab. f 
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for multiple detections of the same source te.iellol l2003). 
If the source shows long-term, erratic variability, it is ex- 
cluded from the list of SN candidates. For the candidates 
which passed this selection and that were centered on the 
host galaxy nucleus (~ 30%) we found that 40% were ac- 
tually SNe while the remaining ones were still AGNs. To 
take this into account, in all the calculations nuclear can- 
didates were given a statistical weight of 0.4. Note that 
with the progressing of the monitoring, AGN contamina- 
tion will continue to be reduced. 

In the 5 fields discussed in this paper, we have de- 
tected 40 candidates. These are listed in Tab l2l where 
we report the candidate designation (col. 1), co- 
ordinates (cols. 2-3), and apparent V magnitude 
at discovery (col. 4), the apparent V magnitude of 



the host galaxy (col. 5), the offset (in arcsec) from 
the host nucleus (col. 6), the host photometric red- 
shift (col. 7), the number of individual exposures 
in which the object has been seen (col. 8) and the 
epoch of first detection (col. 9). In the last col- 
umn we report the classification of the candidate. 
Candidates located in the host galaxy nuclear re- 
gions are labelled SNAGN and are given a smaller 
weight as described above. 

For 9 of the candidate we obtained spectroscopic ob- 
servations which were used to derive spectral classifica- 
tion and redshift (Tab. |3J). Eight of these objects were 
confirmed as SNe, classified as type la (3), type II (4) and 
type Ic (1). One object was found to be a Seyfert 1 galaxy. 
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Table 3. SN candidates with spectroscopic classification 



Designation 


Type 


Redshift 


Field 


Reference 


SN1999ey 


Iln 


0.094 


J1888 


IAUC 


7310 


SN1999gt 


la 


0.275 


AXAF 


IAUC 


7346 


SN1999gu 


II 


0.149 


AXAF 


IAUC 


7346 


SN2000fc 


la 


0.430 


J1888 


IAUC 


7537 


SN2001bd 


II-L 


0.096 


Field2 


IAUC 


7615 


SN2001gh 


II-P 


0.159 


J1888 


IAUC 


7762 


SN2001ip 


la 


0.536 


AXAF 


IAUC 


7780 


SN2002cl 


Ic 


0.072 


13Z3 


IAUC 


7785 


13Z3-K 


Syl 


0.362 


13Z3 







3. The Galaxy Sample 

The estimate of SN rates in the local Universe relies on 
the charac terization of the galaxy sample which have been 
searched ijCappellaro et all Il997|) . For each galaxy one 
needs to know the distance, which enters the computa- 
tion of the surveillance time, and the integrated luminos- 
ity, which is used as a normalization factor. Indeed, it 
has been demonstrated that the SN rate scales with the 
size of the parent stellar po pulation as measured from the 
integrated blue luminosity l|Cappellaro et al~lll9 931. For 
nearby galaxies the relevant information can be readily 
retrieved from published catalogues, while this is not the 
case when one goes to larger distances. 

In an attempt to follow the same approach for our 
intermediate redshift SN project, we exploited the B, V, R 
images obtained during the search to measure magnitudes 
and colors of the galaxies detected in our fields and we 
used them to derive distances and absolute luminosities 
through the SEP fittin g photometric redshift technique 
llRolzonella et al. lEoOrl 

With this aim we selected from the image archive, for 
each field and band, the exposures obtained under the best 
conditions, in particular those with seeing < 1", which 
were stacked together using the SWARP package by E. 
Bcrtin (cf. http://terapix.iap.fr/). This produces a sensi- 
ble gain in the S/N ratio: while the limiting magnitude 
of a typical search image is V ~ 24.5, the same for the 
stacked image is V ~ 26 (3 — <r point source). 

From the V images, which compared to the other 
bands benefit from longer cumulative exposure times, 
we built the galaxy catalogue in cluding all sources with 
SExtractor stellarity index < 0.9 l|Arnouts et alJl200l[) . A 
further selection is performed excluding galaxies fainter 
than R = 21.8. While this limit was originally chosen to 
improve our confidence in the photometric rcdshifts (see 
next), it has also the advantage that up to this magnitude 
the S/N is good enough to guarantee a clean star-galaxy 
separation and a good photometric accuracy. For these 
selected galaxies, B and R magnitudes were eventually 
measured adopting the same aperture as defined in the V 
image. 

P hotometric redshifts wer e estimated using the hyper-z 
code l)Bolzonella et al. Il2000l) . This program searches for 



the best match between the measured colors of galaxies 
and the values in a grid created from a library of spec- 
tral energy distribution (SED) templates for different red- 
shifts. It has been shown that with a proper sampling of 
the SED, even with broad band filters, galaxy redshifts 
can be measured with a n RMS error as small as a ~ 0.05 
( Bolzonclla et al. 2000). In our case, however, with obser- 
vations in only three bands we cannot match this level of 
accuracy. In particular, we note that one of the most sig- 
nificant features in the galaxy SED, the 4000 A break, falls 
redward of the R band for redshifts larger than z > 0.8. 
This means that the redshifts derived for more distant 
galaxies are very uncertain. 

On the other side we notice that, because of the observ- 
ing strategy and limiting magnitude, the SNe discovered 
in our search are all at z < 0.8 with a peak of the dis- 
tribution at z ~ 0.3. To remove as much as possible the 
contamination of distant galaxies erroneously estimated at 
low redshift, we removed from the sample galaxies fainter 
than R = 21.8. This roughly c orresponds to the magni- 
tude M* at a redshift z = 0.8 ijWolf et al. Il2003h. where 
M» i s a parameter of the Schechter function ijSchechterl 
Il97fih which is used to fit the galaxy luminosity function. 
With this choice, galaxies with redshift > 0.8 contribute 
about 20% of the total sample luminosity. 

At the same time, low luminosity galaxies at low red- 
shift are removed from the sample. However at z ~ 0.3, 
the average redshift of our search, they contribute to only 
1/4 of the total luminosity and they are expected to give 
a small contribution to the SN productivity. Indeed we 
found that out of the 40 SN candidates, only 2 (AXAF-H 
and AFAX-A, which were not included in the computa- 
tions) were discovered in galaxies fainter than R = 21.8. 

The final galaxy catalog, which counts about 11300 
galaxies, was cross-checked with NED in order to assess 
the accuracy of the photometric redshifts. We found that 
spectroscopic redshifts were available for 324 galaxies (in- 
cluding 118 galaxies of the field J1888 for which spec- 
troscopic redshift were kindly provided by P.-A. Due, in 
advance of publication) . With these data we built the his- 
togram of Az = (z p h — 2 sp )/(l + z sp ), the differences be- 
tween spectroscopic and photometric redshifts, which is 
shown in Fig. ^ (dashed line). It turns out that the aver- 
age difference is (Az) = 0.01 with a RMS error a — 0.10, 
which is both consistent with our limited SED sampling 
and sufficient for our statistical analysis. An independent 
check on the accuracy of our photometric redshift was 
made possible by the fact that one of our fields (AXAF) 
partially overlaps with the Chand ra Deep Field Sout h cov- 
ered by the COMBO-17 survey l|Wolf et al. 112003^ . The 
distribution of Az for 1375 galaxies which we have in com- 
mon is also shown in Fig^ (thick line). The average dif- 
ference (Az) = 0.00 with a RMS error a — 0.11 is very 
similar to those measured in the comparison with spectro- 
scopic redshifts. 

The redshift distribution of our galaxy sample is shown 
in Fig. El As expected, given the adopted limiting mag- 
nitude, the number counts peak at about z ~ 0.3. For 
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Fig. 1. Distribution of Az = (z p h — z sp )/(l + z sp ), the dif- 
ference between our estimate of the photometric redshift 
and a) spectroscopic redshifts, for the 324 galaxies of our 
sample with known spectroscopic r edshift (thin line) and 
b) COMB017 photometric redshifts lWolf et al. I l|2003ft for 
the 1375 galaxies which are in common with our sample 
(thick line) 
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We need to estimate the time during which a SN hosted 
in a given galaxy remains brighter than the search detec- 
tion threshold which is usually referred to as control time. 
This depends on i) the search threshold and detection ef- 
ficiency; ii) the SN absolute magnitude and luminosity 
evolution, Hi) the distance. 

There are several concurrent factors determining the 
detection efficiency. For a given instrument and fixed ex- 
posure time, the key factors are the observing conditions, 
in particular sky transparency and seeing. There is also 
some dependence on the characteristics of the particular 
sky field, due to the disturbing presence of bright stars 
and nearby galaxies. The position of the candidate within 
the parent galaxy plays also a role. 

While the sky transparency of a specific observation 
is measured by the photometric constants, the impact of 
other factors was evaluated performing a number of arti- 
ficial star experiments. In these simulations we placed a 
number of synthetic stars with a given magnitude m, and 
the PSF deduced from field stars, distributing them in 
different galaxies of the field. The position of the artificial 
stars within the host galaxy was chosen randomly, assum- 
ing a gaussian distribution centered on the galaxy nucleus 
and a FWHM equal to that of the host galaxy. Then, 
the synthetic frames were processed through our search 
pipeline and the detection efficiency e(m) was computed 
as the ratio between recovered and injected stars. 

In Fig. |3| we plot an example of detection efficiency 
function for one of the observation of the field AXAF. 
As a result of these numerical experiments, which will be 
described in detail in a forthcoming paper (Riello et al. 
2004), we found that the most critical parameter is the 
seeing. In particular we found that, taking as reference the 
magnitude at which the detection efficiency is 50%, this is 
^1.2 mag fainter when the seeing is 0'.'65 compared with 
a seeing of 1'.'3. 

The other factor entering the computation of the con- 
trol time for a given epoch is the SN luminosity evolution 
and its absolute luminosity. The prescription we used in 
our calculations is described in the following. The appar- 
ent V light curve for each SN was computed taking into 
account the luminosity distance (estimated from the pho- 
tometric redshift of the host galaxy), the time dilation and 
the i^-correction according to the following relation: 



Fig. 2. Redshift distribution of our galaxy sample (line, 
left-side scale) and SNe host galaxies (dots, right-side 
scale). For the latter, candidates coincident with the host 
galaxy nucleus have been counted with a 0.4 weight. 

comparison, in Fig. we have also plotted the redshift 
distribution of the SN candidate host galaxies (dots) with 
their statistical errors. 

4. Detection Efficiencies and Control Time 

The computation of rates requires the definition of the 
time interval during which the events could be detected. 



(1) 



where, for the i-th galaxy at redshift z,*, to = t/(l + Zi) 
is the galaxy rest frame time, M§ N (io) is the SN light 
curve in B absolute magnitude, fJ,(zi) is the galaxy dis- 
tance modulus 3 , i^£y(io) is the B to V if-correction and 
Ay is the galactic extinction in V (the upperscript SN 
indicates SN type dependent quantit ies). 

Following lCappellaro et all (|l997l ). we considered four 
basic SN types, namely la, IIP, IIL and Ib/c. Guided by 
the general theoretical interpretation, type IIP, IIL, and 

3 Hereafter, we assume a standard flat cosmology Om = 0.3 
and fi A = 0.7 and h = H /75 
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Fig. 3. Example of SN detection efficiency curve, plotted 
as a function of apparent magnitude for one of the epochs 
of field AXAF (seeing 0"9). 



Ib/c will be collectively called later core collapse (CC) 
SNe. B ab solute magnitude s and light curves were also 



taken as in lCannellaro et al The iv^y-correction 

for each SN type at different phases was calculated us- 
ing a sample of spectra with good S/N ratio available 
in the Asiago-Padova SN archive. Spectra were redshifted 
by the proper amount and synthetic photometry was com- 
puted using the synphot package in IRAF. The choice of 
the Kbv correction is motivated by the fact that, for the 
average redshift of our search, (z) ~ 0.3, this minimizes 
both the uncertainties and the time dependence of the 
^-correction. Note that for the same reason we use as 
reference B light curves of local SNe. 

We introduce here the effective control time defined, 
for each observation of a given galaxy, as the control time 
weighted by the galaxy luminosity in units of 10 10 Lbq: 



5. SN Rates 

For a given galaxy and SN type, the rates are derived 
by dividing the number of observed events by the total 
control time: 



„SN 



(1+*) 



GT 



■SN 



(3) 



where the factor 1 + z corrects the rate to the rest 
frame. 

The average redshift of the SN search (zsn) (he. of 
the galaxy sample) is given by the average of the galaxy 
redshifts weighted by the effective control time: 



A'* 



E z i CT 



;SN 



(zsn) 



jY 



(4) 



where N is the number of galaxies of the sample. It results 
that (zi a ) = 0.32, (z cc ) = 0.26 for la and CC SNe respec- 
tively. The lower (z) for CC SNe is obviously due to the 
fact that, on average, these are intrinsically fainter than 
la. 

Considering that we have a sample of galaxies with 
a wide spread in redshift, to derive an estimate of the 
SN rates we have computed the expected SN counts as a 
function of redshift, derived by summing the contribution 
of individual SN types and using different assumptions for 
the rate evolution with redshift. For the latter we assume 
a power law dependence rgN(^) == ?sn(1 + where 7g N 
represent the local (z = 0) rate for a given SN type and 
a is the evolution index (a = indicates no evolution). 
Using this model, the SN rate evolution parameters can 
be derived from the best fit between the expected and 
observed redshift distributions. 

In practice, we derive the total expected SN detection 
N exp (z) as the sum of the expected number of la and CC 
events: 



N^(z) = N°: P (z) + N^(z) 



(5) 



The expected redshift distributions of type la and CC SNe 
are given by: 



■SN 



U 



SN 



(m) e(m) dm 



(2) 



where rf N (m) is the time a SN in the i-galaxy stays at 
a magnitude between m and m+dm, and can be computed 
from relation QJ, e(m) is the detection efficiency at the 
given magnitude and Li is the luminosity of the i-th galaxy 
in units of 10 10 blue-band solar luminosities. 

Finally, for the given galaxy the total control time 

SN 

of the search campaign, CT, ; , is obtained by properly 
combining the control time of individual observations (cf. 
ICappellaro et al.lll999j) . 



i=l 



(6) 



(7) 



where the sums are extended over the n galaxies in a given 
redshift bin z and the effective control time for core col- 
lapse SNe is computed as follows: 



CTT(z) = / Ib/C CT? /C (z) + f u CT?(z) 



(8) 



where we introduced the relative rates of type II, fu = 
r n/ r cc an d type Ib/c, /ib/ c = 1 ~ fu- Here we make the 
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Fig. 4. Galaxy lu minosity density ad different re d shifts. 
Data are from lIMadau. Pozzetti. fe Dickinsonl Il998t 
Wolf et al. | 



. J2003L Blanton et al. I l2003t iNorberg et al. I 
2002^ . The line represents the best fit using a (1 + z) a 



power law (see the text for details). 



further assumption that the relative rates of different CC 
types do not change wit h redshift and are equ al to the 
local ones, i.e. fa = 0.8. llCapoellaro et alJ ll999h . 

As a first step, we have tested the null hypothesis 
that the SN rates per unit blue luminosity do not change 
with redshift (a/ = a cc = 0) and they are equal to 
the local values, namely r ° = 0.20 ± 0.05 \h 2 SN u] 4 and 
r° c = 0.47 ±0.19 [fr 2 SNu] (ICappellaro et al.lll999h . Under 
this hypothesis, in our galaxy sample one would expect 
13.2 events, which is significantly smaller than the 28.2 
actually observed. This already points out that the global 
SN rate must increase rapidly with redshift. 

In our case, since we do not have the spectroscopic clas- 
sification for all the SN candidates, we cannot directly use 
equation to derive individual SN rates. Unfortunately, 
with the limited statistics available to date and the in- 
complete spectroscopic classification, we cannot determine 
rates and evolutions of both SN la and CC SNe from our 
data alone and we need to use other data sets. For this 
purpose we have exploited the publi shed measurements 
of SN la rate at diffe r ent redshifts felanc et al. ll20o4 
Madgwick et al. I l20o3 lHardin et al. I l2000t IPain et al.1 
2002HTonrv et al. 1120031) to fix the evolution of SN la and 
to derive the CC rate. 

One further complication is that some of the high red- 
shift SN la measurement s were given only p e r unit volume 
llMadewick et al. I l2003t IPain et al. 1 120021: iTonrv et al.l 
2003|. To convert these numbers into rate per unit lu- 
minosity (SNu scale), we need consistent estimates of the 
luminosity density at different redshift. For this purpose 
we collected from the literature recent measurements of 
the B luminosity density for the redshift range relevant for 
our SN search (Fig.0J. As can be seen, the measurements 
show a relatively high dispersion, particularly large at high 




0.2 0.4 0.6 
redshift 

Fig. 5. Measurements of the SN la rate at differ- 
ent r^d^hifts^^el^ences are as follows: filled triangle 

- llCappellaro et all Il999l) . filled upsi de-down tria ngle 

- felanc et al. I l2004f) , filled square - iMadewick et al~1 
120031). empty square - l|Hardin et al. I 20001). filled circl e 



l|Tonrv et ai~ll20q3h. filled diamond - JPain et al. 1120021) . 



empty diamond - l|Gal-Yam. Maoz. fe SharonlfcOO? ). The 
line is a fit with the power law in (1 + z) (see the text for 
details). Because of the large errorbars the Gal- Yam et al. 
estimates have not be used for the fit. 

redshift. The best fit with a power law p = po(l + z) ap 
gives po = 1.18 x 10 s [/iL BQ Mpc~ 3 ] and a p = 1.0. This 
has been used, in particular , to c onvert the rate per unit 
volume of IPain et al~l l)2002|) and lTonrv et al. 1 1)2003|) . 

The published values of the SN la rate (in SNu) are 
plotted in Fig. [31 which convincingly shows a redshift evo- 
lution. This can be fitted with a power law relation with 
< = 0.18 ± 0.04 h 2 SNu and a Ia = 1.5 ± 0.6. 

Assuming that the SN la rate evolution is known, we 
can use the observed SN counts to constrain the evolution 
of core collapse rates. To this end, it is convenient to re- 
write eq. Q as follows: 



(i + zY 



(1 + (*cc» a < 



-CT, (z) 



(9) 



4 1 SNu = 1 supernova century 1 10 10 L 



B0- 



which is referring the rate to the average redshift 
(z C c) = 0.26 in our case of the galaxy sample. 

By means of a maximum-likelihood method we search 

(z ) 

for the values of rcc and a cc which give the best fit 
between observed and expected distributions as a function 
of redshift. Resulting confidence levels are shown in Fig. O 
We found that r cc {z = 0.26) = 1.45±^;|g h 2 SNu and a cc = 
2.9^2 g where the quoted errors are the 1-er confidence 
level. Clearly the parameter a cc describing the CC-rate 
evolution is poorly confined by our data alone, but the fair 
agreement with the measurement in the local Universe (cf. 
Fig. |HJ is contorting. 

The comparison between the observed and expected 
distribution for the best fit parameters is shown in Fig.0 
We notice that in our search about 2/3 of the observed 
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Fig. 7. Observed (points with statistical errorbars) and 
expected (solid line) redshift distribution of SN counts. 
The short-dashed and long-dashed lines are the expected 
type la and CC SN counts respectively. 

SNe are expected to be CC and only 1/3 SN la. This is 
consistent with what we have found when spectroscopic 
confirmation was available (cf. Table OJ- 

The conclusion is that the CC SN rate at (z) ~ 0.3 
appears to be a factor 3 higher than in the local Universe. 

6. Uncertainties 

The errors quoted before for the rate and evolutionary 
index parameter are purely statistical. Even though they 
are large, systematics errors need also to be considered. 
Indeed, although we made our best to exploit the signifi- 
cant observational effort, this first estimate of the CC rate 
at z > should be regarded as preliminary in many re- 
spects. In the following we highlight what we believe are 
the most severe caveats. Note that we will not discuss 
here systematics which originate from the uncertainties 
on the SN properties (absolute magnitudes, light curves, 
intrinsic dispersion). For them, in fact, we have made the 
same assumptions as for the computation of local rates 



l|Cappellaro et alJll999^l and therefore they are expected 
to cancel out in the comparison. Also we do not address 
the uncertainties in the cosmological model which was 
adopted and we rather focus on other aspects. 

As we mentioned above, the most severe concern is 
the lack of spectroscopic classification for all candidates. 
In particular, we had to account in a statistical manner for 
the expected contamination by AGNs, assigning a weight 
of 0.4 to the candidates coinciding with the host galaxy nu- 
cleus . In order to evaluate the impact of this assumption, 
we have computed the best fit for two extreme cases: a) all 
nuclear candidates are AGNs or b) 80% of the nuclear can- 
didates are indeed SNe. The two extreme cases encompass 
the ranges: ric cc> = 1.3 -r- 1.9 h 2 SNu and a cc = 3.4 2.3 
for the two fit parameters respectively. Therefore although 
important, this uncertainty is not likely to affect our con- 
clusions. 

Another concern is related to the limited accuracy in 
the characterization of the galaxy sample through photo- 
metric redshifts, due to the poor sampling of the galaxy 
SEDs. To check the possible influence of this uncertainty, 
for the 1375 galaxies of the AXAF field which we have 
in common with the COMB017 surveys we computed 
the expected number of events alternatively using our 
o wn estimates f or the galaxy redshifts and those reported 
in IWolf et al~l l|2003h . We found that to match the SN 
number count in the AXAF field, using as reference the 
COMB017 redshift estimates the rate rc C has to be 
higher by ~ 10% compared with our own redshift esti- 
mates. This is a small difference which is completely hid- 
den by the statistical errors. 

Of more concern is that in the current work we did not 
attempt to correct for the bias due to extinction in the host 
galaxies. This applies to all other estimates of SN rates at 
high redshift published so far. It is well known that in local 
SN searches there is a severe bias for SN detection in spiral 
galax ies which are not observed face-on ((Cappellaro et alJ 
1999). This is attributed to the concentration of dust in 
the plane of the spiral galaxy disk which causes a higher 
average extinction for the inclined ones. Because of the 
wavelength dependence of extinction, it is expected that 
the blue photographic surveys which were used to derive 
the estimates of the local rates are heavily affected. For the 
same reason, it is sometimes claimed that CCD searches in 
nearby galaxies do not need to be corrected. This cannot 
be entirely true because the bias is seen even in the list of 
events di scovered in the last few y ears mainly by CCD SN 
searches (|Cappellaro et al. 1120041) . But even if we accept 
this claim, we stress that if high redshift galaxies have 
the same dust content as local ones then, just because 
of redshift, host galaxy extinction should cause a similar 
bias for blue band searches in nearby galaxies as for red 
searches in high redshift galaxies. Unfortunately, current 
data on high redshift SN searches do not allow to measure 
the size of this bias. 

At the moment, we can only add that any correction 
is likely to increase, possibly even significantly, the CC 
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Fig. 8. CC SN rate with redshift. The dot is 
the estimate of the l ocal CC SN rate from 
ICappellaro et al.l (|l999^ whereas the square is 
the new measurement derived in this paper. The 

solid line shows the deduced evolution with redshift. The 
dashed line shows the deduced type la rate evolution (see 
the text for details) . The thick solid line is the CC SN ra te 
evolution predicted by model M2 of ISadat et al. 1 1 1998f) . 



estimate reported here. In summary none of the reported 
uncertainties, even if important, seems to undermine the 
main conclusion reached here, that a significant evolution 
of the CC rate must be present even for this small look 
back time. 



7. Discussion 

When we compare our result with the models we find 
that, in the redshift range we are exploring, the cur- 
rent predictions indicate a much shallower evolution 
of the CC rate llSadat et al. I Il998t iMadauet al. I Il998l 



ICalura fc Matteucci 1 120031 iKobavashi et al. I l200C|) than 
actually observed.. 

As an ex a mple, in FigJS] we show the predictions of 
ISadat et al. I l)l998|) in the most favorable case (model 
M2 which corresponds to a higher SFR at high redshift). 
Although at a 2-ct level the model is consistent with our 
measurement, we remark that taken to face value the ob- 
served rate at z=0.26 is a factor 2 higher than the model. 
This becomes even more significant if we consider that, 
due to the lack of extinction correction, our estimate is 
expected to be a lower limit to the actual CC rate. 

On the other hand, the CC rate evolution is directly 
related to the adopted history of SFR for which there are 
many new estimates. Indeed, considering the short evolu- 
tionary time scale of the CC SNe progenitors (< 5 x 10 7 
yr) and under the assumption that the initial mass func- 
tion (IMF) and the mass range of the CC progenitors do 
not change significantly in the redshift range of interest, 



there is a simple direct relation between the SFR and the 
CC rate, namely: 



J^I ct>(M)dM 



^ u M<p(M)dM 



(10) 



where ip is the SFR at the given redshift, 4>{M) is 
the IMF, M L - M v is the mass range of the IMF and 
Mf c — M^ c the mass range for the CC SN progenitors. In 
particular, adopting a Salpeter IMF in the range Ml = 0.1 
to M v = 125 M Q and Mf c = 8, M£ G = 50 M Q for the 
lower and upper limits of the mass of CC SN progenitors 
we derive r cc ~ 0.007^- 

Usually the cosmic SFR is reported per unit of co- 
moving volume and hence, for the comparison with SN 
rate evolution, it is convenient to translate the latter to 
the same unit using the known evolution of the luminos- 
ity density (Fig. 0J. With this conversion, the CC SN 



1+0.8 
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X 



rate per unit volume results r^ c (z = 0.26) = 2.2 
10 -4 /i 3 yr _1 Mpc -3 (whereas the local value translates in 
r^(z = 0) = 5.5 ± 2.2 x 10" 5 h^T^Mpc' 3 ). Finally, 
using eq. 10, we converted the CC rate in SFR, obtain- 

-3 



3.1±i o x 10~ 2 h 3 M 



yr -1 Mpc 3 (and 



/i 3 M Q yr" 1 Mpc^ 3 for the local 



ing tp(z = 0.26) 
ij>(z = 0) = 7.9±3.1xl0 
value). 

Recently, iFuiita et al. ] 1120031) have published an esti- 
mate of the SFR based on the Ha luminosity density at 
z ~ 0.24. This is shown in Fig|5] where are also reported 
othe r estimates of the S FR at different redshifts (adapted 
fro mlFui ita et al." (I2003T) , see references therein) . As noted 
bv lFuiita et al. I l|2003h there is a systematic difference in 
the SFR deduced from the Ha luminosity density (filled 
symbols) compared with that derived from the UV lumi- 
nosity density non corrected for extinction (empty sym- 
bols), with the former being significantly smaller at all 
redshift. However, it has been shown that this disagree- 
ment can be removed with a proper extinction corr ection 
l|Rosa-Gonzalez et al. ll2002t iHinnelein et al. Il2003h . 

As seen in Fig. 9, our measurement i s in e xcellent 
agreement with the value of IFuiita et al. I l)2003|) and in 
general with the rapid SFR evolution deduced from Ha 
luminosity density. Our conclusion, from the evolution of 
the core collapse SN rate is that at at a redshift z = 0.26, 
that is at a look back time of 2.8 Gyr with the adopted cos- 
mology, the SFR per unit comoving volume was three time 
higher than in local Universe. We stress again that, be- 
cause our measurement of the CC SN rate is not 
corrected for extinction, this is likely to be a lower 
limit (cf Sect. 6). 5 



5 After submission of this article, two preprints have 
been posted reporting new estimates of the SN rate 
evolution up to redshift 1.6, as par t of the G reat 
Observatories Origins D eep Survey dStrolger et alTI 



120041 ; iDahlen et al. 1 1200 J) . The main caveat is that, 
like in our case, they could obtain spectroscopic clas- 
sification only for a fraction of the candidates. In any 
case, it turns out that their estimate of the CC SN rate 
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Fig. 9. We compare our estimate of the SFR at redshift 
z = 0.26 (big fille d dot) with the recent estimate of 



iFuiita et al. I l|2003|) based on theHa luminosity density 
at z = 0.24 (filled square). Also shown are estimates of 
the SFR at other redshifts based either on measurements 
of the Ha (fill ed symbols) or of the UV luminosity density 
(adapted from lFuiita et al"l l|2003f) and reference therein). 
Also plotted is the value derived from the local CC rate 
(also big filled dot). 
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